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Abstract 



The purpose of this paper is to find and analyze plane symmetric, 
static(non static) solutions in Hofava- Lifshitz gravity. We discussed two 
versions of Horava gravity. First we showed that if the detailed balance 
principle have considered, there are both static and non-static solutions. 
We show that in static case there are two family of solvable models which 
either of them has a well defined EOS, in analogous to the perfect fluid 
' solutions in GR. In non-static case we find a family of solutions. Some 

physical properties of these solutions was discussed. Secondly we investi- 
gated the plane symmetric solutions for a new modified version of Hof avaa 
gravity [18], which has the new terms inserted action in it. 



1 Introduction 



The Einstein's equations for hypersurface-homogeneous space-time are reduced 
to a system of ordinary differential equations. At least for the spatially-homogeneous 
' case, they form a well-posed Cauchy problem [T] and, although they have not 

been completely integrated, their qualitative properties have been discussed in 
many papers: (see Wainwright and Ellis [5] for an extensive survey of the re- 
sults, and e.g. [3] Ryan and Shepley , MacCallum [3], Bogoyavlenskii [5], and 
Rosquist and Jantzen [6] for useful earlier reviews). Methods of dynamical sys- 
tems theory which proved fruitful in elucidating these properties lead to ways 
of restricting the general case to more readily solvable subcases and thence to 
new exact solutions (see e.g. Uggla et al. [7] for a summary). Nearly all of 
these methods were initially developed for use in the spatially-homogeneous 
case (cosmologies, for brevity), and we shall therefore describe the methods in 
this context although they can be adapted to the G3 on T 3 and H 3 on V3 cases 
also (e.g. for an ortho normal tetrad method for G3 on T3 see Harness The 
number of freedom degrees , i.e. the number of essential arbitrary constants re- 
quired in a general cosmology for each Bianchi type, has been studied by Siklos 
[5] and MacCallum 10 , Wainwright and Ellis [H]). The residual set of ordinary 
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differential equations which are going to be solved can be formulated in various 
ways [2J, Wainwright and Ellis [2J. One can use a time-independent basis as 
in and parameterize the components g a p in suitable ways. This is called the 
metric approach. One may then choose spatial coordinations and a new time 
coordination r. 

The main alternative to the metric approach is the ortho normal tetrad ap- 
proach using a tetrad basis. The two ways are closely related when variables 
are chosen as just described above [TT], in both approaches a suitable choice 
of lapse (or, to include the G3 on T 3 case, slicing gauge), or directly of inde- 
pendent variable, may decouple and simplify the equations; for details see e.g. 
Jantzen [12], Uggla et al. [7]. A power-law lapse, i.e. a product of powers of 
the dependent variables, an idea introduced in Bonanos [T3], may be useful, as 
may an intrinsic slicing (e.g. making a product of invariantly-dcfincd metric 
components the independent variable). 

Recently, a power-counting renormalizable, ultra-violet (UV) complete theory 
of gravity was proposed by Hof ava in [IH [15j [16j EZ] ■ Although presenting an 
infrared (IR) fixed point, namely General Relativity, in the UV the theory pos- 
sesses a fixed point with an anisotropic, Lifshitz scaling between time and space 
of the form x % — > I x % , t — > £ z t, where £, z, x % and t are the scaling factor, 
dynamical critical exponent, spatial coordination and temporal coordination, 
respectively. But recently Bias, et.al listed some inconsistencies of the Hofava- 
Lifshitz gravity as a complete description of Quantum gravity. They discussed 
the consistency of Hofavas proposal for a theory of quantum gravity from the 
low-energy viewpoint. They uncover the additional scalar freedom degree aris- 
ing from the explicit breaks the general covariance and its followed properties. 
Their analysis was done both in the original formulation of the theory and in 
the Stuckelberg frame. A unusual feature of the new mode is that it satisfies 
a first order (in time derivatives) equation of motion . In linear aproximation 
this extra freedom degree manifested only around non-static spatially inhomo- 
geneous backgrounds. They found two serious problems associated with this 
mode. " First, the mode develops very fast exponential instabilities at short dis- 
tances. Second, it becomes strongly coupled at an extremely low cutoff scale" . 
They also discussed a version of Hofavas gravity with projectable condition and 
stated that this version is a certain limit of the ghost condensate model. The 
theory is still problematic since the additional field generically forms caustics 
and, again, has a very low strong coupling scale. Also they clarify some sub- 
tleties that arise from the application of the Stuckelberg formalism to Hofavas 
model due to its non-relativistic nature. 

Due to these features, there has been a large effort on examining and extend- 
ing the properties of the theory itself [HI [201 EU [221 H31 [23 HH [2H [301 
I3"T1 [231 [231 152") . Additionally, application of Hof ava-Lifshitz gravity as a cos- 
mological framework gives rise to Hofava-Lifshitz cosmology, which is going 
to lead from interesting behavior [33j [3H [35j [36j [37]. In particular, one can 
examine specific solution subclasses [551 I3H HO]) the perturbation spectrum 
[UJ [H gl gl i51 SSI E7], the gravitational wave production gSJ 05], the mat- 
ter bounce [SOI EQ [S2J, the black hole properties [51 (Ml 023 US EZl GSH GE3 , 
the cosmic string solutions [5D] the dark energy phenomenology EJJ E21 E3 El] > 
the astrophysical phenomenology [65j [66] etc. However, despite this extended 
research, there are still many ambiguities if Hofava-Lifshitz gravity is reliable 
and capable from a successful description for the gravitational background of 
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our world, there still many ambiguities. [2"31 [Ml 125] . 

In the present paper we would like to obtain plane-symmetric perfect fluid so- 
lutions in the framework of Hofava-Lifshitz gravity. At first we review such 
solutions in both static and non-static cases in the framework of general relativ- 
ity. After a brief review of Hofava-Lifshitz gravity, we obtain and analyze plane 
symmetric, static Taub-like solution. We extend our study to the Purely time 
dependent solution. Also we discussed some simple solutions for a new version 
of modified Horava theory which was recently proposed by Bias , et al. 

2 Plane-symmetric perfect fluid solutions in GR 

2.1 Static solutions 

The plane-symmetric static perfect fluids with a prescribed equation of state 
fi = /x(p) are given by [57J HH] 

ds 2 = -~e 2v dt 2 + z 2 (dx 2 + dy 2 ) + zFiz^dz 2 (1) 



where 



2V z 3 „ , 

= I - k qP — = -22i/' (2) 



fJ-{p)+P F 

F' = -K 0l i(p)z 2 (3) 

For a given function fi = fi{p), the differential equations determine F(z) and 
p = p(z) and from p(z) then v = v(z). Equations (2,3) lead to the condition [69 



z 2 + G G + z 3 
G z 3 - G 

kqp(z) 

So one may prescribe G(z) and then obtain p and v — v[z) as linear integrals. 
For an equation of state p = (7 — l)fi the function G was found by Collins [70] . 
A solution for p = /i/3 was obtained by Teixeira et al. [71]. The solution for 
/i = const was given by Taub [67] and by Horsky [72] in terms of hypergeometric 
functions. Some other special solutions have been given by Davidson [73] . Static 
plane-symmetric perfect fluids occur also as subcases of the static cylindrically- 
symmetric solutions. 




2.2 Non-static solutions 

Besides the classes of solutions which were described in a, several other classes 
have been found by making special assumptions for the metric or the equation 
of state. For p = /z, Tabensky and Taub 68 reduced the field equation to a 
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single linear differential equation, 

ds 2 = r 1,2 e n {dz 2 - dt 2 ) + t{dx 2 + dy 2 ),t > (4) 

+«*>*+ *^ w 

<r,tt + t V,t - ct zz = (6) 

Ko p = K /i = * 1/2 e-V 2 t - a 2 ,), (7) 
Tabensky and Taub [55] also gave the special solution 

cr = a log(t) + /3 arccos(-) (8) 

ft = 2(a 2 + /3 2 ) log(i) + 2/3 2 log(l - Z —) + 4af3 arccos(-) (9) 



This solution represents the asymptotic Robertson- Walker Space-time. These 
functions are particular solutions which are inhomogeneous and anisotropic but 
toward to an homogeneous spacetime for large times. The solution depends 
on two parameters a and (3. For a 2 + (3 2 — §. We get Robertson- Walker 
space-time asymptotically. The variables z and t are restricted to z 2 < t 2 . 

For further solutions see Collins and Lang [74] , Goode [75] , Carot and Sintes 
[76] . Bray [77], Tariq and Tupper [78], Shikin [79], Gupta and Sharma [80]. We 
want to obtain the analogous of the metric solution (3) in Horava theory. We 
use from this gauge, since firstly it has a well known simple GR analogous and 
secondly since the field equations has only one gauge function and also this 
solution is spatially Ricci flat, i.e R = and several terms in Horava action will 
be vanished, especially the Cotton tensor Cy = 0. 

3 Review of Horava-Lifshitz gravity with detailed 
condition 

Following from the ADM decomposition of the metric [ST], and the Einstein 
equations, the fundamental objects of interest are the fields N(t,x),Ni(t,x), (t, x) 
corresponding to the lapse , shift and spatial metric of the ADM decomposition. 
In the (3 + l)-dimensional ADM formalism, where the metric can be written as 

ds 2 = -N 2 dt 2 + g l3 (dx l + N l dt)(dx j + N j dt) (10) 

and for a spacelike hypersurface with a fixed time, its extrinsic curvature K%j is 

Kij - 2^(fl« - ViiVj - VjNi) 

where " dot" denotes a derivative from " t" and covariant derivatives defined with 
respect to the spatial metric gij , the action of Hof ava-Lifshitz theory for z — 3 
is 

S= I dtd 3 Xy^N{C K - C v ) (11) 

J M 
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we define the space-covariant derivative on a covector Vi as VjU; = diVj ~ TyVi 
where T\. is the spatial ChristofFel symbol. g is the determinant of the 3-metric 
and N = N(t) is a dimensionless homogeneous gauge field. The kinetic term is 

C K = \o K = -^ji A', l\' J - XK 2 ) 

Here Ni is a gauge field with scaling dimension [Ni] = z — 1. 
The potential term ZV of the (3 + l)-dimensional theory is determined by the 
principle of detailed balance [3J, requiring £y to follow, in a precise way, from 
the gradient flow generated by a 3-dimensional action W g . This principle was 
applied to gravity with the result that the number of possible terms in Cy are 
drastically reduced with respect to the broad choice available in an potential is 

C v = aedjCV - a^R^jR™ 1 + <> ; /,',., /'"' - A ^ X _\) R2] + ^ ~ 3Aw) (U) 
The coupling constants cti define by 

Ct4,A w 



Q'2 



3A- 1 

OLA = 

8 

K 2 fl 

2v z 



Where in it Cy is the Cotton tensor [3] which is defined as, 

T kB? 



Following 38] we can write the action as 



S = 



J dtdx 3 (£ +£i) (13) 



C Q = ^N(l( Kij KV - XK 2 ) + K ^ 2{ 8 ^ R ^ Al) ) (14) 



4 Exact solutions 

We shall restrict ourselves to situations where the space time has plane symme- 
try.the problem is invariant under transformations of the form 

x — > x + a 

y y + b 

x — > x cos(9) + y s'm(6) 
y —> ycos(8) — xs'm(6) 
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The metric of a spacetime which admits plane symmetry may be written as [81] 

ds 2 = _ e 2F dt 2 + l( e 2H (da ,2 + dy 2, + tf) (16) 

Where F, G and H are just functions of z and t alone. By substituting (16) in 
(14,15) finally we obtain the following form of action 

S = fdtJ d 3 x[^e 2H+G+F [^e- 2F [(l \)G 2 + 2(1 - 2A)ff 2 - 4XGH] + g( ff 3X) (A^-R - 3A 2 J (17) 
+ K 32( 1 (1 _~ 3 ^ ) fl2 " ^y^e" 4G ((2i? 2 + H" - H'G' f + (2ff 2 + 2H" - 2H'G') 2 )]] 

Where in it /' = f = * and the spatial Ricci Scalar for metric (15) reads 
as, 

R = 2c 2 e- 2G (3H' 2 + 2H" - 2H'G') 

By varying G, H and f , we obtain the equations of motion. Considering the 
case where 3A — 1 > 0, one can see, from the modified Friedmann equations in 
the context of the Horava theory, a negative A is required. In Ref. [55]) the 
authors point out that a positive A can be obtained by making an analytical 
continuation for parameters. In this work we emphasized on the3A — 1 < 0. 
Choosing A like above avoids us from a negative cosmological constant term 
which is so harmful for the cosmological evidences. 

4.1 Static Taub-like solution 

In this section we want to obtain a Taub like solution (1) as metric functions 
(2,3) with action (17). Comparing two metrics (16) and (l)we observe that we 
must take the metric functions as H = ]og(cz),G = log(cy^^y) and F = v{z). 

In this case the scalar Ricci is R = 2^ and the full action (17) convert to the 
next functional: 



S = J dtd 3 x$(v,F(z),F'(z)) (18) 
<%, F(z),F'(z)) = -±-c 3 e» K 2 n 2 z- 3 - 5 F(z)-°< 5 (4\ - 1)~\3X - l)" 1 ]T f n z n (19) 



n=0 



32 

The v equation of motion is: 

Uz n - (20) 

n=Q 

Where the Functional coefficients /„ were defined by: 

f = (24A 2 - 11A+ l)F 2 (21) 

/j = (-144A 2 + 90A - U)FF' (22) 

h = (24A 2 - 11A + 1)F' 2 (23) 

h = 8A W (4A-1)F (24) 

U = 8A W (4X-1)F' (25) 

/ 6 = 12A,„ 2 (1-4A) (26) 
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Just F(z) field equation is more complicated, we do not write it here. We men- 
tioned here that the general solution for metric function F(z) and the energy 
density function p for action (19) is: 

(27) 

F(z) = -e- 1 /(a, 6; z)[12A(4A - l)I(c, d; z) + (1 - 3A)ci] 

h = 8X(z-3)- z + 7 

I(a,p;z) = { h a z?dz 



K p = e- 1 z- 2 [h a z b (12X(4A - l)J(c, d; z) + (1 - 3A)ci) + 12A(4A - l)h c z d I(a, b; z)} (28) 
Where 

704A 2 - 400A + 55 
°~ ~ (8A- l)(24A-7) 
„ 4A - 1 

b = — 

c= 16 



24A-7 
(4A-l)(8A-3) 



(8 A- l)(24A-7) 
40A- 11 



d = 2 "24A-7 
e= (-1 + 3A)" 1 

the integral can be calculated by using the Hypergeometric functions atH 

h c z d dz = —(7 - 24A) C F[1 + d, -c, 2 + d, \ _ ^ x l ] (29) 

Without overloading the paper with lengthy formulae, the pressure p(z) and 
metric function v(z) can be calculated from (2,3). It was argued by Horava 
that in the infrared limit the higher order terms become negligible and so one 
may expect to recover GR if the parameter A flows to 1 in that limit. But 
Bias et al [TH] argued that this expectation is incorrect: the explicit breaking 
of general covariance leads to the appearance of an extra degree of freedom 
in the infrared which becomes strongly coupled when A approaches 1. More 
precisely, the additional mode is weakly coupled only in a narrow window at low 
energies. This window depends both on A and the parameters of the background 
geometry; it shrinks to zero when A — > 1 or when the background curvature 
vanishes. Thus the solution has no familiar GR limit. Also as stated by Bias et 
al |18) , the only consistent solution for this theory must be time dependent and 
inhomogeneous. Another reason for considering the plane symmetric metrics 
goes to the Mukohyama work [S5]. He considered the version without detailed 
balance condition with the projectablility condition and addressed one aspect 
of the theory: avoidance of caustics for constant time hypersurfaces. He showed 
that there is no caustic with plane symmetry in the absence of matter source if 
A ^ 1. If A = 1 is a stable IR fixed point of the renormalization group flow then 
A is expected to be deviated from 1 near would-be caustics, where the extrinsic 
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curvature increases and high energy corrections become important. Therefore, 
the absence of caustics with A ^ 1 implies that caustics cannot form with this 
symmetry when the source is absent. He argued that inclusion of matter source 
will not change the conclusion. He also argued that caustics with codimension 
higher than one will not form because of repulsive gravity generated by nonlinear 
higher curvature terms. These arguments support Mukohyama conjecture that 
there is no caustic for constant time hypersurfaces. Our exact plane symmetric 
solution endorsement this conjecture directly. 



4.2 Other special solutions 

Case A = 1/3: 

As the observation in [16] the value A = 1/3 corresponds to the action being 
invariant under an anisotropic conformal (Weyl) symmetry. The case A = 1/3, 
corresponds to an anisotropic Weyl-invariant theory, is another relevant value 
for which the equations of motion simplify and hence admit quite explicit so- 
lutions. Recently |84) derived general static spherically symmetric solutions in 
the Hof ava theory of gravity with nonzero shift held. These represent hedgehog 
versions of black holes with radial hair arising from the shift held. For the case 
of the standard de Witt kinetic term (A = 1) there is an infinity of solutions 
that exhibit a deformed version of re paramctrization invariance away from the 
general relativistic limit. Special solutions also arise in the anisotropic confor- 
mal point A = 1/3. This is a simple reason that why we investigated solution 
with A = 1/3. 
The field equation is: 

A w z 3 = ^(zF)' (30) 

The general solution for this differential equation is: 

F(z) = ^zJ + Cl /z (31) 

To find the unknown metric function v we must use (2,3). By solving the simple 
differential equations for v(z),p( t z) and by substituting ji{z) from (3), after a 
simple integration it leads to thi 
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^) = l/8(-^£+ 8 ^)(32) 

Q 

p(z) = -0.3735£A- (33) 

n 

= Jc[^P(-0.25, 1.0308, 0.35354=) - 3.1768ciP(0.75, 1.0308, 0.35354=) (34) 

s/ci Jc{ 



+c 2 Jo[£Q(-0.25, 1.0308, 0.35354=) - 3.1768cic 2 Q(0.75, 1.0308, 0.35354=) 



•n = J^n (-i 2 + 8ci)(c 2 Q(-0.25, 1.0308, 0.35354=) + P(-0.25, 1.0308, 0.3535-^)) (35) 



3 P(a, b, x), Q(a,b,x) are the associated Legendre functions of the first and second kind, 
new coordinate £ is defined as § = (3A m 2 4 + 8C1) 1 / 2 
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The expression of the metric function v is volumeless and we do not write 
it here. We must mention here that this spacetime represents a new exact 
solution for GR with a complicated implicit EOS p = p(fi) which could be 
derived by eliminating the z coordinate between the both pressure and energy 
density expressions given by the above functions. 
Case A = 1/4: 

In this case the differential equation (20) is: 

^z 2 {F'f + zFF' + ^F 2 =0 (36) 

This equation has a simple solution 

F{z) = ^ (37) 

this is the special case of solution (28) with the A w — 0. We can write the 
following expressions for fi(z),p(z): 



p(z) = fx(z)(4: - VT7tanh^-(c 2 - logz)) (39) 
by eliminating the z between (35,36) we can write the EOS as: 



p(p) = fx(4 - VT7tanh^(c' 2 + 1/4 log p)) (40) 

Where c' 2 — c 2 — l/41og(ci//co)- 

The metric function fj,(z) in this case has simple expression as: 

v{z) = log(^/ 2 cosh ^Qog z - 4)) (41) 
Thus the metric (16) transforms to the following form: 

ds 2 = z 2 [-z- 1 ' 2 cosh(-^J(logz - c' 2 ))dt 2 + dx 2 + dy 2 + ^dz 2 ] (42) 



4.3 Purely time dependent solution 

If we consider the non-static, plane symmetric solutions with the metric anstaz(4). 
In synchronous time t, the Plane ADM metric has Ni — , where fi is the 
usual Tabensky and Taub metric function, we assumed that the metric function 
51 = fi(t). Since The metric (3) is Spatial Ricci flat, and further the Spatial Ricci 
tensor and consequently the result Cotton tensor are vanished, then the full 
curvature dependence part of the action (13), i.e the L\ — O.The easiest way 
to obtain the solution for the full Lagrangian is to substitute the metric anstaz 
(3) into the action, and then vary the function Q. This is a valid procedure in 
analogous with cylindrically and spherical kinds. The result reduces Lagrangian 
to an overall scaling constant, is given by 

c = J,r 3/2 e- a 4 + f 5/2 {-f- 5,A + t-^ti? \{\ + mf] z^r x Al) (43) 
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Perhaps surprising power 5/2 can be understood by observing that this form 
of C is invariant under the reflection, we assumed that the metric function 
Q = fi(t). The remaining equation is obtained by variation of Lagrangian (15). 
The field equation is: 

d (A = bV - fe n (44) 



dV dt 
where 



(45) 



ds 2 = " 1 ^^ U v^ ^ -^l± (d z i -di 2 )+ t(dx 2 + dy 2 ) (47) 



(3A-1)(A-1) 
The general solution for this ODE is 

nM , , (ggjj ^ 9) sec 2 (lV8 Cl 6 - 9(c 2 - log(t)) 

= iog( ) ( 46 ) 

Thus the general solution for metric (3) with Lagrangian (15) in context of 
Horava theory is: 

(8ci6 - 9) sec 2 (|V8cib-9(c 2 - log(*)) 

The interpretation is as follows: we choose d as the proper distance from z = 
along x, y, i constant. The density vanishes at z = oo and at z = 0. The plane 
2 = behaves like a hard wall. It can be shown that all time-like geodesies 
oscillate between two z planes determined by initial conditions. No time-like 
geodesic ever touches the wall. We can, imagine the space-time for z < to 
be the vacuum solution. Strictly speaking that one doesn't know the matching 
condition for singular hypersurfaces but the type of divergence from both sides 
of z = is exactly the same. We conclude that the only way to have a static 
solution is to introduce an additional boundary condition.We shall now show 
that in asymptotic limit, space-time is homogeneous, and for some of the param- 
eter values, Robertson- WalkeiQ. All we have to do is to carry out the inversion 
for comoving coordinatations for large ts . The calculation is straight-forward 
but rather long, so we omit it. Therefore as similar to the GR, we get that 
for any group of fluid world-lines, the space-time is homogeneous. If the four- 
velocity is irrotational, it can be written in terms of a scalar function a(z,t) 
as u n — (— a a <j a ), and Modified Einsteins field equations for a stiff fluid read 
p = p. 



5 Plane symmetric solution for modified Horava- 
Lifshitz gravity 

In this section we attempt to find analytic solutions for the modified Horava- 
Lifshitz gravity proposed recently by Bias, et. al [18| . For this we choose the 
action to be 

S = J d 3 xdtN^{a{K lj K lJ - \K 2 ) + pC t] C rj + £R+ a^a.a 1 )) (48) 
° i; = — 

4 Asymptotic Robertson- Walker Space-time 



10 



The same analysis was performed by this action and by imposing static gauges 
for Lapse and the metric functions (with a vanishes shift function) by Kiritsis is 
Spherical symmetry [83]. and we are going to look for plane symmetric solutions 
with zero shift of the form 

ds 2 = ~N{z) 2 dt 2 + dz 2 + u(z) 2 (dx 2 + dy 2 ) (49) 

For such static-gauge the kinetic terms do not contributes the equations of 
motion. We insert this metric ansatz in the action we have: 



d 3 xdtL 
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L = 2f3Nu' 2 ( — - \) 2 + 2i(u' 2 + 2m/') + a lU 2 -— 

5.1 Solution in the absence of the Cotton tensor 

There is only one simple solution in case (3 = 0. In this case the contribution 
from Cotton tensor disappeared. This solution is written as: 

N{z) ={i cot ^ i{ci \!% z) m 

u(z) = c 2 , v ' (51) 

5.2 Constant (spatially) curvature solutions 

We know that constant curvature solutions form are a rich family of investiga- 
tions both in GR and f(R) metric gravity.If we accept that there is a Cosmologi- 
cal constant term in the universe which dominated in epoch, thus all the models 
of the universe imposed that the field equations govern the spacetime manifolds 
must be asymptotically de Sitter. This means that the spacetime posses a con- 
stant Ricci scalar. Also in the context of the string theory the existence of the 
constant negative curvature solution (AdS) is fundamental. Far away from any 
of this reasons as an attempt we can suppose that there is a constant spatial 
Ricci scalar solution for our field equations in this modified form of the Horava 
theory. This family of exact solutions can be obtained by supposing that for 
this solution the Ricci scalar R is constant. Solving the field equation for u(z) 
we obtain: 

p2A I R 

U( - Z) = V"R cosh ^ + V 2A Z ) ( 52 ) 

for finding the lapse function N(z) one must solve the following differential 
equation 

^ + ^ + 2(log(u(z))Y)=g(z) (53) 

V N 

«)- ,{z) 



2aiu{z) 2 

/(*) = 2/?«'(*)(f -2(^1) 2 ) 
4 u(z ) 
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This is a standard form of a Riccati equation. There is not found a general 
solution for it, yet for non vanishing values of g(z). But if we set (3 = then 
this differential equation solves by simple methods. The final result is: 



N(z) = c 2 [2Atanh 2 (W — z) + 4\/2ARci tanh(W — z) + ARc\] (54) 



6 Conclusion 

The Horava-Lifshitz gravity is a power-counting renormalizable theory, which 
has anisotropic scale between space and time in the UV limit, and thus breaks 
the Lorentz invariance.By applying this theory to plane symmetric spacetimes, 
one finds a set of simple differential equations. This system could be solved 
analytically and two family of static Taub like solutions are obtained. In both 
families the cosmological constant term is positive. We mentioned here that we 
do not have any analytical continuation on the parameters of model. We observe 
that in contraction of the cosmological models based on the Horava theory, there 
are no trouble problematic features of model in the special case A = 1/3, 1/4. 
So the negative cosmological constant problem can be solved successfully in the 
plane symmetry. In this case, the system has both a state and a non static one. 
The former corresponds to the two families corresponds to the A = 1/3, 1/4. The 
first solution A = 1/3 is a generalization of the Taub fluid solution in GR.But the 
EOS of the fluid in this case is so complicated and contains Legendre functions. 
Once the equation of state w — £ of a perfect fluid reaches to infinity. If 
we rewrite the field equations in terms of a first order system of autonomous 
equations, this stable critical point coincides with the saddle one. Thus the 
stable state is broken. It may be expected that the equation of state may change 
with the energy scale. Actually, if A = 1/3, this solution represents the role of 
an infinite charge in the Kehagias-Sfetsos (KS) black hole[85]. We notice that if 
g^ is a solution then ^g^ is also a solution, whenever if is constant. From now 
on it shall be understood that all line elements can be multiplied by a constant 
conformal factor. The non static solution recovers the GR solution in IR regime. 
It is homogenous and asymptotically Robertson- Walker Space-time in comoving 
coordinations. In spite of the GR stiff matter solutions, unfortunately there are 
not many non-static solutions with only one Killing vector.We may can state a 
theorem as: 

Theorem: If the metric (4) satisfies the vacuum field equations for a Horava 
type of gravity, then the perfect fluid field equations for a stiff fluid are satisfied 
by the metric (4) which differs from (3) by the some simple (but not so trivial 
as the GR) in metric function Sl(t, z) respectively. 

We lead proof of this theorem to another work. The technique described by the 
above theorem can be applied to many solutions: the class of vacuum metrics 
(4) belongs to those metrics with an orthogonally transitive Abelian groups. 
Moreover, since the metric function g xx = g yy = t 2 does not change under these 
generation techniques, the same functions a and SI can be used for all vacuum 
solutions obtained from the same vacuum seed. A different way of looking at 
the class of solutions covered by this Theorem is for starting from a stiff fluid 
solution, performing a transformation to a vacuum solution, applying a soliton- 
generation technique and then going back to the stiff fluid: one then may speak 
of a solution describing solitons travelling in the background of a (particular) 
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stiff fluid of higher symmetry. Of course, since a and f2 need not be changed, it 
suffices to immediately transform only when the vacuum part of the metric. 
Also we derived two family of exact solution for the modified version of Horava 
theory which recently proposed by Bias et al.lt was shown that there is one 
solution which it's spatial curvature is constant. Another solution is obtained in 
the absence of the Cotton tensor. 
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